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EFFECTS  OF  CAVITATION  ON  UNDERWATER  SHOCK  LOADING  - 
PLANE  PROBLEM,  Part  2 

1.  Introduction 

Development  of  a  finite  element  fluid  code  for  solving 
plane  problems  involving  shock-induced  cavitation  was  described 
in  Ref.  1.  The  program  uses  displacement  potential  as  the 
dependent  variable  as  described  in  Ref.  2.  A  structural  program 
(Ref.  3)  providing  a  plane  model  for  a  ring-stiffened  cylin¬ 
drical  shell  has  recently  been  completed.  In  the  following 
sections,  coupling  of  the  structural  and  fluid  programs  is 
described  and  results  thus  far  are  summarized. 

2 .  Structural  Program 

2.1  Program  Features.  The  structural  program  ( S TRUK 5 ),  was 
developed  by  Jack  T.  Waller,  a  thesis  student  of  mine,  and 
is  described  in  Ref.  3.  A  shell  structure  having  specified 
radius,  area  per  unit  axial  length,  radius  of  gyration,  mater¬ 
ial  density,  and  modulus  of  elasticity  is  modeled.  Trigono¬ 
metric  series  are  used  to  represent  the  radial  and  tangential 
components  of  shell  displacement.  Taking  account  of  symmetry, 
the  half-circumference  of  the  shell  is  divided  into  N  equal 
segments.  The  fluid  mesh  is  matched  with  the  structural 
model,  providing  an  interface  node  at  each  point  of  division. 
Using  fluid  pressures  calculated  at  these  nodes  by  the  fluid 
program,  the  coefficients  for  trigonometric  series  represen¬ 
tation  of  the  interface  pressure  are  found  by  STRUK5. 
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For  each  trigonometric  order  m  there  is  a  term  c^cosmd 
for  interface  pressure,  a  term  amcosm0  for  radial  displace¬ 
ment  and  a  term  bmsinm0  for  tangential  displacement.  Because 
the  structure  is  rotationally  symmetric,  there  is  no  coupling 
between  coefficients  with  different  orders  m.  For  a  given 
m,  differential  equations  governing  am  and  bm  are  coupled. 

Time  integration  is  performed  using  an  explicit  central 
difference  algorithm.  A  chosen  time  step,  which  may  be 
suitable  with  a  particular  fluid  mesh,  may  result  in  instabi¬ 
lity  for  structural  modes  with  large  m.  Accordingly  the 
program  limits  the  number  of  structural  modes  retained.  For 
this  purpose  the  extensional  modes  are  considered  separately 
from  the  flexural  modes. 

2 . 2  Nonstructural  Mass.  The  shell  structure  of  a  submarine 
typically  has  a  mass  no  greater  than  20  percent  of  the  dis¬ 
placed  water.  The  remaining  mass  is  connected  to  the  shell 
through  relatively  flexible  structure  and/or  isolation  mountings. 
Although  the  early-time  response  to  shock  loading  is  little 
affected  by  this  additional  internal  mass,  its  omission  from 
the  structural  model  will  be  significant  at  later  times. 

Since  cavitation  effects  on  shock  response  may  occur  at  late 
times,  it  was  deemed  important  to  include  the  additional  mass. 

A  rigid  body  displacement  of  the  shell  can  be  described 
by  trigonometric  terms  of  order  m=l.  For  modeling  purposes 
the  spring  forces  between  shell  and  contents  are  distributed 
in- corresponding  fashion.  Thus  there  is  a  set  of  radial 
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forces  varying  as  cos8  and  a  set  of  tangential  forces  varying 
as  sinQ.  This  produces  equations  of  motion  in  which  the 
displacement  of  the  nonstructural  mass  couples  only  with 
the  shell  displacements  described  by  a^  and  b^. 

Modification  of  the  equations  of  motion  to  include  this 
nonstructural  mass  ("contents”  mass)  was  accomplished  in 
the  following  way.  The  notation  of  Waller  (Ref.  3)  is  used 
with  the  following  additions . 

uc  =  contents  displacement  (in  8  =  180°  direction) 

m  =  contents  mass 
c 

kc  =  common  stiffness  of  radial  and  tangential  spring  sets 
The  increments  to  elastic  potential  energy  and  kinetic  energy, 
AU  and  AT,  are  given  by 


AU 

*  1i*c[<a1-uc>2  ♦  Cbj-u^2] 

(1) 

AT 

= 

(2) 

When  these  additions  are  included  the  motion  equations  for 
a^  and  b^  each  contain  both  a^  and  b^.  Algebraic  manipulation 
allows  the  equations  to  be  written  in  a  form  suitable  for 

explicit  time  integration  as 

2  .  c, 

(l+2z)a1+  ^■(l+z)(a1+b1)+^^[(l+z)a1-zb1-(l+2z)uc]=(l+z)^  (3) 
R 

2  z  c 

(l+22)b1+27(l+z)(a1+b1)+^||s-[-za1+(l+z,)b1+(l+2z)uc]  =  — ~  (4) 

R 

The  additional  equation  of  motion  for  the  contents  mass  is 
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0 


(5 


mu  +  k  ( 2u  -a.,  +b. )  = 
c  c  c  c  1  1 


The  contents  mass  mc  is  determined  within  the  program  from 
the  stipulation  that  shell  plus  contents  be  neutrally  buoyant. 


Instead  of  specifying  kc  directly,  the  fixed-base  frequency 
f  of  the  contents  is  specified  and  kc  is  found  from 

k  =  f  2/<8ir2m  ) 
c  c  c 


(6) 


2.3  Curvature  Correction.  When  the  just-described  process 
for  including  nonstructural  mass  in  the  model  was  first 
implemented  it  was  discovered  that  effects  of  shell  rotatory 
inertia  appeared  in  the  equation  for  rigid  body  motion  of 
the  shell  (obtained  by  putting  b^=-a^  in  the  m=l  equations). 
Since  this  behavior  is  clearly  incorrect,  the  equations  were 
carefully  re-examined  to  discover  the  cause.  Two  related 
errors  were  discovered.  Both  appear  first  in  my  own  notes 
and  are  clearly  my  responsibility,  not  Waller's. 

The  first  mistake  is  that  the  expression  used  for 
rotational  velocity  u  of  the  shell  is  incomplete.  The 
correct  form  is 


.  1  3w  <r 

“  ‘  K  TRF  "  K 


(7) 


The  second  term  is  omitted  in  Ref.  3.  The  other  mistake  is 
in  the  expression  for  shell  bending  curvature  k.  The  correct 
form  is 
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In  Ref.  3  the  second  term  in  parentheses  is  +w.  The  latter 
form  is  correct  only  if  the  ring  deformation  is  inextensional . 

i 

The  reason  that  these  errors  were  not  discovered  during 
initial  testing  of  Waller's  program  is  that  they  have  negli¬ 
gible  effects  on  results.  Nevertheless,  the  corresponding 
changes  have  been  made  in  the  structural  program.  The  new 
version  of  the  program,  which  also  includes  the  nonstructural 
mass  as  detailed  in  Section  2.2,  is  called  STRUK6. 

Revisions  to  the  equations  of  Ref.  3  resulting  from  the 
corrections  to  rotational  velocity  (Eq.  7)  and  curvature 
(Eq.  8)  are  given  in  Appendix  A. 

3 .  Fluid  Program 

Principal  features  of  the  fluid  program  were  described 
in  Ref.  1.  Additional  details  are  given  in  Ref.  4,  including 
discussion  of  the  staggered  solution  procedure  which  is  used 
to  couple  fluid  and  structure  for  a  time  marching  solution. 

4.  Results 

Two  kinds  of  results  have  been  obtained  from  program 
tests  conducted  to  date.  The  first  kind  give  no  new  infor¬ 
mation  concerning  cavitation  effects,  but  simply  demonstrate 
that  the  coupled  response  of  the  structure  is  in  agreement 
with  known  analytic  results.  These  data,  which  afford 
supporting  evidence  concerning  the  validity  of  the  modeling, 
are  summarized  in  Appendix  B. 
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Results  of  the  second  kind  concern  effects  of  cavitation 


upon  the  severity  of  shock  induced  structural  loads.  For 
the  cases  studied  to  date  results  are  summarized  in  Ref.  4 
and  quoted  below. 


Examination  of  results  from  a  large  num¬ 
ber  of  cases  studied  up  to  this  time  reveals 
that  the  ratio  of  shock  pressure  increment  to 
hydrostatic  pressure  (pg/p^)  must  exceed  three 
times  the  ratio  of  decay  length  to  structure 
radius  (L/R)  in  order  to  induce  cavitation. 
Cavitation  begins  about  one  decay  length  in 
front  of  the  structure  at  about  the  time  when 
the  initial  scattered  wavefront  reaches  this 
location . 

Extreme  values  of  bending  moments  and 
circumferential  compressive  force  in  the  ring 
have  been  found  to  occur  before  the  shock 
front  has  completely  traversed  the  structure. 
Since  cavitation  has  not  yet  begun  during 
this  interval,  the  extreme  structural  loads 
are  not  affected  by  cavitation.  There  is,  to 
be  sure,  a  shock  pulse  generated  by  cavity 
closure.  The  initial  magnitude  of  this  pulse 
is  of  the  same  order  as  hydrostatic  pressure. 
After  travelling  to  the  structure  the  pulse 
is  significantly  attenuated  and  the  resulting 
loads  are  not  significant. 


5.  Conclusions 

Cases  investigated  thus  far  have  shown  no  instance  in 
which  the  occurrence  of  cavitation  increases  the  maximum 
stresses  in  the  structure.  A  wider  range  of  encounter 
parameters  needs  to  be  investigated  before  it  is  possible 
to  assert  with  confidence  that  cavitation  induced  by  structural 


response  to  underwater  shock  will  not  increase  the  likeli¬ 
hood  of  structural  damage  or  collapse. 

Considering  the  related  phenomenon  of  bulk  cavitation, 
evidence  to  date  does  not  rule  out  the  possibility  that  it 
may  intensify  shock-induced  structural  loading.  The  finite 
element  method  should  be  used  to  investigate  effects  of 
bulk  cavitation  on  structural  loads . 
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APPENDIX  A 


Corrections  for  Equations  of  Reference  3 

Corrected  forms  for  equations  of  Ref.  3  based  on  expressions 
for  u)  and  k  given  in  Eqs .  7  and  8  of  this  report  appear  below. 
Equation  numbers  are  those  of  Ref.  3. 


5U  .  ^Xn2(nan  +  bn)(n5an  +  5bn)  +  ^(an+  nbn)(5an  +  n6bn)  (10) 

Rs 

dT  *  ^p(ARd9)(w2  +  v2)  +  |pIRde(jJ|-  -  0 


<5T  *  irpAR£(an<San  + 

(1  +  n2z)a„  +  nz’b* 
n  n 

nzan  +  (1  +  z)bn  + 

n*  *  1  *  (n2  -  l)2z 
n2 


6A}  +  +  6n)(n5*n  +  SbJ 


+  (0(1  +  n*z)an  +  (0  n(l  +  n2z)bn 


(02n(l  +  n2z)an  +  (0*n2(l  +  z)bn 


3  +  (0  (fl2  -  1)2Z3  =  ~~r 


=  0 


(12) 


(14) 


(17) 


/ 
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Equations  14  each  contain  both  an  and  bR  and  are  thus  unsuit¬ 
able  for  use  with  an  explicit  time  integration  algorithm. 

Some  algebraic  manipulation  produces  the  following  forms 
which  are  used  in  STRUK6 . 


(1  *  2  +  n'2)3n  *  ©’[>  ♦  *  <"’  -  l)!z]a„  *  (£)'(!  *  z)nt>n  .  (1  ♦  z)fn_ 

(1  ♦  z  *  n!z)bn  ♦  (1  -  z  *  2n!z)nan  t  (|)’(1  t  n!z)n!bn  .  -nz^jj- 


APPENDIX  B 


Modal  Frequencies  of  a  Submerged  Cylinder 

It  is  well  known  that  late  time  response  motions  of 
submerged  structures  can  be  accurately  forecast  by  treating 
the  surrounding  fluid  as  incompressible  (Ref.  5).  This 
permits  calculation  of  the  frequencies  for  modes  involving 
either  rigid  body  translation  or  shell  flexure.  Comparison 
between  such  calculated  results  and  observed  frequencies 
of  the  shock  response  motion  of  the  finite  element  model 
provides  a  partial  check  on  the  accuracy  of  the  modeling  and 
the  implementing  programs . 

Figure  1  shows  the  rigid  body  displacement  a-^  of  the 
shell  and  the  displacement  uc  of  the  contents  mass  as  a 
function  of  time  in  response  to  a  shock  with  step  pressure 
rise  of  4MPa  and  decay  time  of  6ms.  In  view  of  the  rapid 
shock  pulse  decay,  the  motion  following  the  first  half 
cycle  may  be  considered  free  vibration.  The  succeeding 
cycle  of  motion  has  a  period  of  77ms  or  a  frequency  of  13.0  Hz. 
For  the  shell  translation  mode  the  added  fluid  mass  is  equal 
to  the  displaced  mass.  Taking  this  into  account,  the  fixed 
base  frequency  f  =  10.0  Hz  for  the  contents  mass  gives  a 
corrected  value  of  13.1  Hz  in  good  agreement  with  the 
observed  result. 

Figure  2  shows  results  for  a  separate  computer  solution 
with  f  =  5.0  Hz.  In  this  instance  the  observed  period  is 
152ms,  corresponding  to  a  frequency  of  6.6  Hz.  This  compares 
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closely  with  a  calculated  value  of  6.5  Hz. 

For  the  flexural  racdes  corresponding  to  the  trigonometric 
orders  n  =  2,  3,  and  4,  the  corresponding  added  mass,  based  on 
maximum  radial  displacement  amplitude,  is  (displaced  mass)/n. 
Table  1  compares  theoretic  frequencies  with  finite  element 
method  (FEM)  results.  It  is  apparent  that  agreement  is 
satisfactory . 


Table  1.  Theoretic  and  Observed  Modal  Frequencies 


Order  n 

Frequency  (Hz) 

Theory  |  FEM 
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